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ABSTRACT 


A distinct visual signature occurs in black holes that are surrounded by optically thin and geomet- 
rically thick emission regions. This signature is a sharp-edged dip in brightness that is coincident 
with the black-hole shadow, which is the projection of the black hole's unstable-photon region on the 
observer's sky. We highlight two key mechanisms responsible for producing the sharp-edged dip: i) the 
reduction of intensity observed in rays that intersect the unstable-photon region, and thus the perfectly 
absorbing event horizon, versus rays that do not (blocking), and ii) the increase of intensity observed 
in rays that travel along extended, horizon-circling paths near the boundary of the unstable-photon 
region (path-lengthening). We demonstrate that the black-hole shadow is a distinct phenomenon from 
the photon ring, and that models exist in which the former may be observed, but not the latter. Ad- 
ditionally, we show that the black-hole shadow and its associated visual signature differ from the more 
model-dependent brightness depressions associated with thin-disk models, because for geometrically 
thick and optically thin emission regions, the blocking and path-lengthening effects are quite general. 
Consequentially, the black-hole shadow is a robust and fairly model-independent observable for ac- 
creting black holes that are in the deep sub-Eddington regime, such as low-luminosity active galactic 


nuclei (LLAGN). 
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1. INTRODUCTION 
1.1. The black-hole shadow 


In the recreational practice of coin rubbing, a sheet 
of paper is placed over a coin, and a pencil or crayon is 
rubbed on the paper. An impression of the coin appears 
in the rubbing. As the trick is repeated, the rubbings 
will look different each time, depending on such factors 
as the type and color of crayon, and the vigor of the 
rubbing; but, as long as a few basic rules are observed, 
the impression of the coin will always be the same. 

An analogy can be drawn with the appearance of a 
black hole, which itself provides a constant spacetime 
structure, but which can be illuminated by external 
sources of electromagnetic radiation, which may have 
a variety of shapes, emit a variety of colors, and be- 
have in a time-dependent manner. Like the coin under- 
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neath a sheet of paper, the black hole itself cannot be 
seen, but its imprint on the light distribution can be 
seen. Recently, the Event Horizon Telescope collabo- 
ration published the first image of a black hole (Event 
Horizon Telescope Collaboration 2019), which exhibits 
this effect of a constant spacetime structure illuminated 
by a time-varying emission region, prompting us to bet- 
ter understand and explain what we see, and what we 
do not see. 

The visual appearance of black holes was first studied 
by Cunningham & Bardeen (1973) and Bardeen (1974) 
for the case of a star orbiting a black hole, as well as var- 
ious related scenarios. This early work revealed the ba- 
sic structure of a thin disk with a gravitationally lensed 
inner region. The first computer-calculated, but hand- 
draw, visualization of a black hole surrounded by a lu- 
minous accretion disk was presented by Luminet (1979). 
These early models focused on geometrically thin, op- 
tically thick accretion disks (Shakura & Sunyaev 1973; 
Novikov & Thorne 1973; Page & Thorne 1974) that have 
so far dominated our perception of black holes in science 
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and popular culture. Such accretion disks may be found 
in black-hole binary systems and quasi-stellar objects 
(QSO’s) (see, e.g., Agol & Krolik (1999)). Those sources 
were found and extensively discussed in the 1960’s and 
1970’s. However, the vast majority of black holes in the 
universe — the silent majority (Falcke 2001) — operate 
at highly sub-Eddington accretion rates. These Low- 
Luminosity AGN (LLAGN) (Ho 2008) are characterized 
by advection-dominated flows that are radiatively inef- 
ficient (see Narayan et al. (1998) for a review), coupled 
to ubiquitous radio jets (Nagar et al. 2001), both of 
which are expected to be geometrically thick and op- 
tically thin at some emission frequency. In contrast 
to Quasars and high-state X-ray binaries, which show 
a plethora of appearances and states, Sub-Eddington 
black holes are reasonably well-behaved (Fender et al. 
2004), and the scaling of their spectral energy distri- 
bution is, to first order, described by a “fundamental 
plane relation” (Merloni et al. 2003; Falcke et al. 2004; 
Plotkin et al. 2012) that represents the ratio between 
optically thick (typically radio) and optically thin (typ- 
ically X-ray) emission as a function of black hole mass 
and accretion rate. The fundamental plane relation also 
determines the frequency at which the emission region 
surrounding a black hole becomes optically thin. This 
model explains why supermassive black holes are opti- 
cally thin somewhere between the radio and NIR fre- 
quency regimes, while X-ray binaries typically become 
optically thin only in the NIR-X-ray regime. For a low- 
power supermassive black hole, like Sgr A* in the cen- 
ter of the Milky Way, optically thin emission occurs at 
mm-wavelengths, and approaches the scale of the event 
horizon (Falcke et al. 1993; Narayan et al. 1995; Falcke 
et al. 1998; Ozel et al. 2000). 

The difference between optically thin and optically 
thick emission regions is conceptually important, be- 
cause optical-depth effects can significantly influence 
the appearance of emission regions around black holes 
(Bronzwaer et al. 2020). The presence of a geometrically 
thick, optically thin emission region directly surround- 
ing a black hole’s event horizon gave rise to the concept 
of the black-hole shadow (Falcke et al. 2000a) that is 
observed when imaging such accreting black holes, and 
whose features depend only on spacetime/observer ge- 
ometry (Johannsen 2016; Johannsen & Psaltis 2010a; 
Bronzwaer et al. 2020). The concept was derived from 
the work by Bardeen (1973), who calculated the optical 
appearance of a black hole in front of a plane-parallel, 
uniformly radiating disk. While this scenario, e.g., a 
black hole in front of a star, is impossible to observe 
with current technology, the key visual element of this 
appearance is the so-called photon sphere - the bound- 


ary of the region around the black hole in which closed, 
circular photon orbits exist. It is this surface that de- 
termines the appearance of the black-hole shadow!. 

Although a black hole’s surrounding emission region, 
which renders the shadow visible in the first place, de- 
pends on plasma astrophysics and can have a highly 
time-dependent appearance (Falcke et al. 2001; Broder- 
ick & Loeb 2006; Noble et al. 2007; Moscibrodzka et al. 
2009), the photon sphere - and thus the shadow itself - 
is as constant as the outline of a coin. As long as the 
emission region is optically thin and geometrically thick, 
a strong visual signature (a dip in intensity) that is coin- 
cident with the shadow should be observable. Moreover, 
the shadow’s radius is proportional to the black hole’s 
mass and, under favorable circumstances, the shadow’s 
shape is informative about the black hole’s spin (see, 
e.g., Johannsen & Psaltis 2010b; Younsi et al. 2016). 
For these reasons, observing the black-hole shadow was a 
key driver for the Event Horizon Telescope (EHT, Event 
Horizon Telescope Collaboration et al. 2019a) and one 
of its clearest predictions (Falcke et al. 2000a; Broderick 
& Loeb 2006; Kamruddin & Dexter 2013; Mosgcibrodzka 
et al. 2016). 


1.2. Criticism of the black-hole shadow 


In the scrutiny of the EHT result subsequent to its 
publishing, criticism has been raised regarding the claim 
of having observed the black-hole shadow. The use 
of the term shadow is also occasionally questioned. 
Bardeen (1973) called his dark region the “apparent 
shape of the black hole". Cash (2002) used the term 
silhouette, in the context of a proposed X-ray inter- 
ferometer to observe accretion disks of X-ray binaries. 
Gralla (2020) used the term “critical curve" to describe 
the projection of the photon sphere's boundary on the 


1 Some of the terms that are used frequently by the Event Horizon 
Telescope (EHT), such as ‘photon ring’ and ‘black-hole shadow’, 
have been assigned multiple interpretations in the literature. In 
order to be consistent with our fellow EHT authors, we will 
henceforth adopt the terminology used by the EHT Gravitational 
Physics Inputs Working Group (GPIWG). In particular, we take 
‘black hole shadow’ (BHS or shadow) to refer to the projection, 
on an observer’s sky, of the region around a black hole in which 
circular orbits exist; this region contains the set of observer view- 
ing rays that cross the black hole’s event horizon. The bounding 
curve of the shadow has also been called the ‘n = co photon 
ring’ and the ‘critical curve’. Our definition of ‘photon ring’ is 
the system of lensed images of the black hole’s environment that, 
to the observer, appear just outside of the BHS. Finally, by ‘cen- 
tral brightness depression’ (CBD) we refer to the characteristic 
central dip in brightness observed in black holes surrounded by 
optically thin accretion disks. As we shall show in this paper, 
systems exist in which the CBD does not match the BHS, but in 
which a distinct visual signature is nonetheless observable at the 
BHS. 


observer’s sky. The term shadow was introduced by Fal- 
cke et al. (2000a) and independently by de Vries (2000) 
a few weeks later, and has been adopted by the Event 
Horizon Telescope (EHT) and a majority of the papers 
that have been published on the topic in the past two 
decades, both on the observational side and on the the- 
oretical side (see, e.g., Grenzebach et al. (2014); Taka- 
hashi (2004); Cunha & Herdeiro (2018)). While precise 
terminology is important, words derived from everyday 
human experience do not have a clear mathematical def- 
inition, nor do they have an equivalent in the realm of 
curved spacetimes. But they do carry a deeper meaning. 
Clearly, there is no real bang in the “big bang”, and it 
is not so clear where the “hole” is in a black hole’; still, 
the general public, as well as the scientific community, 
intuitively know what these terms mean. Hence, we here 
continue to use the term shadow for the geometric shape 
of the projection of the photon sphere on an observer’s 
sky, a shape that corresponds to a characteristic CBD 
in low-power black holes. 

Besides the semantics of our terminology, Gralla et al. 
(2019) pointed out that certain emission models give rise 
to CBD’s whose shapes are determined by the charac- 
teristics of the emission region, rather than the space- 
time geometry; in some of these models, no strong vi- 
sual signature is seen at the BHS. More generally, a 
concern exists that the EHT has produced images of 
the detailed astrophysical circumstances of the accretion 
disk/jet system, rather than of the black-hole shadow 
as an observable that is directly informative about the 
black hole’s spacetime geometry. Additionally, it is 
claimed in Gralla (2020) that the black-hole shadow, 
as an observable feature, is “not generic”, but that the 
photon ring - a term which has been defined in several 
ways, but which we here take to mean: the complex 
system of lensed images of the emission region which 
nearly coincides with the boundary of the shadow - is 
generic. Finally, it is claimed in the same paper that 
“the black-hole shadow does not occur in the EHT mod- 
els of M87*”. Taken together, these claims pull into 
doubt the validity the black-hole shadow, both as an as- 
trophysical observable and as a theoretical concept. It 
is therefore important that they be addressed. 

Gralla et al. (2019) presents several examples of emis- 
sion regions whose CBD’s deviate from the black-hole 
shadow. These may be roughly divided into two cate- 
gories: the first category describes ‘backlit’ black holes 
which are illuminated by distant screens, as in (Bardeen 


? For more on the etymology of the word “black hole", see Herdeiro 


& Lemos (2018). 
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1974), and the second category are optically and geomet- 
rically thin disks. The former class of models cannot re- 
produce the M87* image, while the latter category can, 
provided the disk is viewed in a nearly face-on manner. 
'These models help us to understand why the geometrical 
thickness is a key ingredient for observing the shadow. 
'This will be explored further in the coming sections. 

It is important to note that a wide range of reasonable 
astrophysical models has been investigated by the EHT 
Collaboration (Event Horizon Telescope Collaboration 
et al. 2019b), as well as in (Bronzwaer et al. 2020) — an 
endeavor that we do not seek to repeat here. These sim- 
ulations include jets and disc structures of various sizes 
for a wider range of parameters, and they display a wide 
range of observed source morphologies. Some models, 
observed at a frequency at which the model is (partially) 
optically thick, show no visible shadow at all. Some 
show a partially obscured shadow, others an apparently 
exaggerated shadow. An interesting class of models ap- 
pears to show a dip in brightness that is attributable to a 
direct image of the event horizon (Chael et al. 2021). In 
short, the overall appearance of the source is affected by 
a number of factors related to the emission region. The 
immutable shadow itself, however, can still be observed 
clearly for compact LLAGN sources, at frequencies at 
which the emission region is optically thin. 

In this work, we try to make more intuitively clear 
why, and under which conditions, the black-hole shadow 
is a robust and meaningful feature. We explore the phys- 
ical processes that give rise to the observability of the 
shadow and show why it is that the shadow is in fact the 
more generic feature for geometrically extended emission 
regions. We also investigate the nature of the shadow in 
intermediate cases, such as toroidal and thin-disk emit- 
ters, which cover a range of astrophysical cases. 


1.3. Outline 


In Section 2, we analyze the visual appearance of dark 
objects surrounded by optically thin emission regions 
theoretically, first in a Newtonian approximation (with 
no light bending), and then in a general-relativistic con- 
text; we identify the key physical effects that cause us to 
observe a darkened region on the sky. Section 3 presents 
numerically computed image maps of black holes sur- 
rounded by various types of emitters, which demonstrate 
how the shadow behaves. We summarize our findings in 
Section 4. 


2. RADIATIVE TRANSFER: THEORY 


In this section, we outline a geometrical-optics treat- 
ment of the appearance of an opaque, absorbing, spher- 
ical object located within a luminous, optically thin 
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cloud. We first examine the classical, Newtonian case, 
from which the effect of gravitational lensing is ab- 
sent, and then move to the general-relativistic treat- 
ment. This allows us to better understand which effects 
produce the shadow, as well as which effects are specif- 
ically caused by the bending of light rays in GR. 


2.1. Newtonian case 


Consider a luminous cloud of roughly spherical shape, 
viewed by a distant observer, so that the viewing rays 
are practically parallel upon arrival. A spherical, opaque 
object is added at the cloud’s center. Given a particular 
observer’s viewing direction, a region now exists from 
which no radiation can reach that observer (see Fig. 1). 
In everyday life such a situation, albeit with a more 
complex geometry, occurs if we look at the wick inside 
the flame of a burning candle or a kettle inside a raging 
campfire. It will be seen as a dark contour. 


Figure 1. An illustration of the blocking effect with Newtonian 
physics. No radiation emitted behind the opaque sphere (solid 
white circle) can reach the observer, so that rays that cross the 
opaque sphere (grey and blue lines) have reduced path lengths 
inside an emission region (dashed white circle) than rays that do 
not cross the opaque sphere (white and red lines). 


In Newtonian physics, the observer’s viewing rays 
travel along straight lines, and an observer image is con- 
structed by solving the radiative transfer equation those 
rays: 

a = jv — av, (1) 
where j is the emission coefficient (which is constant 
in our simple model), J, is the ray's specific intensity, 
and o; is the absorption coefficient (which is zero in our 
model). By considering two rays, one of which passes 
through the sphere and the other doesn't (the blue and 


red lines in Fig. 1, respectively), we define the ‘darkening 
factor' x as the ratio of the intensities of the two rays 
as they reach the observer: 


= ie jv ds + Ly sph 

f? jv ds 
where s parametrizes the integration distance along the 
rays, [spn is the specific intensity emitted by the sphere, 
and the position markers a, b, c, d are indicated in Fig. 1. 
Note that x, is frequency-dependent. Therefore, we as- 
sume that the cloud is observed at a frequency v at 
which it is optically thin. 

In the case of a perfectly black sphere, Isp, 7 0 for 
all frequencies. Of course, this dark sphere occurs due 
to infinite redshift at the event horizon in GR, and thus 
would not occur in a Newtonian context. The point of 
this model is only to illustrate the principles that cause 
us to see a shadow. The sphere’s location is near the 
cloud’s center; this is a reasonable assumption, given the 
fact that the sphere, and not the optically thin cloud, 
will dominate the local gravitational field. In fact, the 
dominant central gravitational potential will force ev- 
ery reasonable emission region to either rotate or fall 
rapidly, which will smooth out plasma inhomogeneities 
over time and force the emission region to be symmet- 
ric around the central potential after a few dynamical 
timescales. Since we assume that the emission coeffi- 
cient j of our model is isotropic and constant along the 
rays, fti ds = 1/2 i ds, so that we have 


(2) 


V 


1 

Xbh ~ 2" (3) 
Note that this is an upper bound on xy, in realistic 
scenarios; for very small emitters, grazing rays will have 
much longer path lengths inside the emission region ver- 
sus plunging rays, the arguments mentioned above will 
break, and the darkening will in fact be more severe, as 
the ratio of Eq. 2 diminishes in magnitude. 

Crucially, in the case of an opaque sphere, the dark- 
ening factor will apply if and only if a ray intersects 
the sphere. The darkening factor thus appears as a step 
function in the image, in other words, the darkened re- 
gion has a ‘sharp edge’. 

To sum up, in the Newtonian case, we should expect 
an optically thin cloud of arbitrary shape, with a black, 
opaque sphere near its center, to display a “shadow”, a 
sharp-edged, darkened — but not dark — circle whose flux 
density is roughly half of the flux density from the sur- 
rounding regions. The reason we see this darkened circle 
is because the sphere blocks radiation from an approx- 
imately cylindrical region from reaching the observer. 
We shall call this the ‘blocking effect’. 


We note that the above considerations hold for emit- 
ters of many different shapes (e.g., toroidal), as long 
as the sphere is completely occluded by the emitter; 
only the darkening factor is then affected by the in- 
stantaneous, detailed state of the emission region, while 
the outline of the black-hole shadow remains constant. 
Consider, for example, a toroidal emitter with a radius 
equal to the radius of the central sphere, viewed from 
within the equatorial plane. In this case the shadow 
will be visible, however, if we now move the observer, 
so that the torus is viewed in a face-on manner, the 
blocking effect vanishes and we do not see a shadow, 
but the hole in the torus. Given a significantly flat- 
tened emitter, the blocking effect is conditional. In the 
extreme case of a purely two-dimensional emitter, such 
as the Novikov-Thorne disk (Shakura & Sunyaev 1973; 
Novikov & Thorne 1973), the blocking effect is also sup- 
pressed, and the outline of the sphere could at most be 
partially revealed in the Newtonian case (see Section 
3.2). 


2.2. General-relativistic case 


How does the preceding discussion change if we move 
to a general-relativistic treatment? In Einstein’s gen- 
eral theory of relativity (GR), the structure of spacetime 
is expressed using the metric metric tensor g,,. Light 
rays move along null geodesics; in a powerful gravita- 
tional field (such as that near the center of the Kerr 
metric considered in this paper), null geodesics deviate 
from straight lines, and become curved. A light ray's 
direction of travel is characterized by the wave vector 
k", which is null-normalized: k^k, = 0. By parallel- 
transporting k” along itself, we obtain the geodesic 
equation: 

k^ Vk" — 0, (4) 
where V,, is the covariant derivative. 

In order to create a virtual image of a black hole, 4 
is solved for initial (camera) rays, and then, given a 
radiative model for an emitting source (e.g., an accre- 
tion disk or jet), radiative-transfer calculations are per- 
formed along the rays. To do so, we employ Lorentz- 
invariant versions of the intensity Z = [,,/v°, the invari- 
ant emission coefficient J = j,/v?, and the invariant 
absorption coefficient A = a,v. We may then solve the 
invariant radiative-transfer equation along a ray: 


dL 

oo = 5 

D7 — AL, (5) 
and simply transform back to the usual quantities after 


integration. 

A key change that is introduced with respect to the 
Newtonian case is the gravitational lens (i.e., the bend- 
ing of light rays by gravitational fields), which, in the 
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case of a black hole, exaggerates the size of the shadow. 
This can be seen by examining Fig. 2: the black hole’s 
gravitational field focuses the rays close to the unstable- 
photon region, so that two rays that were widely sepa- 
rated at the observer may both plunge. 

Another key consequence of the switch from a New- 
tonian to a relativistic picture is the introduction of the 
event horizon, the one-way causal membrane that sepa- 
rates the black hole’s interior from its exterior, meaning 
that no information can be transmitted from the for- 
mer to the latter. Although the horizon (R = 2GM/c? 
in the Schwarzschild case) is a black sphere, it is actu- 
ally the unstable-photon region (at R = 3GM/c? in the 
Schwarzschild case), which perfectly captures all rays 
emanating from a distant observer that cross its bound- 
ary. It plays the role of the dark sphere we introduced in 
the Newtonian scenario. Crucially, that spheroidal re- 
gion appears to us as sharp-edged; the reason for this is 
that whether or not a ray crosses this region is a yes/no 
question. Rays that intersect the unstable-photon re- 
gion (and only those rays) are blocked from traveling 
further, thus reducing their path lengths, particularly 
through the emission region of a typical accretion disk 
(see the blue and red rays in Fig. 2). Note that, as the 
black hole is engulfed by its emission region, the ob- 
served dip in intensity will never be complete, i.e., the 
shadow will never be perfectly dark. 

Because the rays are bent toward the singularity by 
the black hole’s gravitational field, rays that ‘graze’ the 
photon sphere begin to ‘loop around’ it; in the extreme 
case of a photon that ends up exactly on the pho- 
ton sphere, it will continue to do so indefinitely. We 
shall call this effect, which occurs for rays that pass 
near the boundary of the unstable-photon region, the 
‘path lengthening effect’. For optically thin sources, 
the path length of a ray through the source is a proxy 
for the intensity picked up by that ray. Thus, the 
path-lengthening effect enhances the source’s brightness 
around the edge of the shadow. We note that the path- 
lengthening effect occurs both inside of and outside of 
the shadow, i.e., for optically thin sources, the boundary 
of the black-hole shadow appears brighter to us than its 
center (note that this effect can also result from emis- 
sivity profiles that increase toward the center; see, e.g., 
Fig. 4 in Narayan et al. (2019)). The blocking effect, 
however, breaks this apparent symmetry, and renders 
the interior of the shadow darker than its exterior. 

Summing up, GR implies two significant changes, 
namely: i) The size of the shadow is increased due to 
gravitational lensing, and ii) The brightness of the cloud 
just outside the shadow (as well as just inside of it) is 
increased due to the path-lengthening effect. 
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Figure 2. A number of null geodesics in the equatorial plane of 
the Schwarzschild spacetime. The red and white geodesics do not 
cross the photon sphere, while the grey geodesics do, causing their 
path length to be reduced, and hence, the observer to see less flux 
in the corresponding viewing directions, as in the Newtonian case 
(Fig. 1). The red rays, which barely escape, are the set of rays with 
at least one turning point. Their paths are lengthened, causing 
the characteristic increase in brightness observed just outside the 
shadow. The solid grey circle marks the unstable-photon region; 
the solid black circle marks the event horizon. 


We emphasize that, while these two effects signif- 
icantly alter the black-hole shadow - increasing the 
shadow’s size and amplifying the source’s brightness 
near its boundary - the primary reason for why we see 
a shadow is still the blocking effect, i.e., the reduced 
intensity from rays that intersect the unstable-photon 
region. That remains true, even if Fig. 3 shows that in 
the case of GR, light can reach the observer from 'be- 
hind’ the black hole due to gravitational lensing. In the 
Newtonian case, the volume behind the absorbing ob- 
ject appears completely dark; in GR, it does not. How- 
ever, the light from ‘behind’ the black hole that reaches 
the observer, will appear to originate from outside the 
shadow to the observer, and only enhance the contrast 
of the shadow further. 

Finally, we note that in optically thick, geometrically 
thin emission regions, such as the Novikov- Thorne thin- 
disk model, neither the blocking effect nor the path- 
lengthening effect can play a role, since, for these mod- 
els, a ray's path length inside the emitter is always zero. 


3. RADIATIVE TRANSFER: NUMERICAL 
CALCULATIONS 


Our next step is to produce images, using simple emis- 
sion models, to investigate the effects described in the 
previous section in a realistic scenario. Since the ap- 


Y [Re] 


Figure 3. Density map of null geodesics that reach a distant 
observer and arrive nearly parallel (below the bottom of the im- 
age), in the Schwarzschild spacetime. Unlike the Newtonian case 
shown in Fig. 1, flux may reach the observer from everywhere 
outside the event horizon (marked by the white dashed circle; the 
black dashed circle represents the unstable-photon region). How- 
ever, a region is still observed from which lower-than-average flux 
reaches the observer. 


pearance of an arbitrary emission region surrounding 
a black hole cannot generally be calculated in closed 
form, we turn to a numerical approximation. We em- 
ploy the general-relativistic radiative transfer (GRRT) 
code RAPTOR (Bronzwaer et al. 2018). In this section, we 
will present image maps of several different source mod- 
els, in order to investigate the behavior of the shadow. 


3.1. Spherical emitters 


Figure 2 shows the geometry of null geodesics for a 
distant observer looking at a Schwarzschild black hole. 
In order for the observer to see anything, an emission 
region must be present that intersects at least some of 
the viewing rays. The simplest possible case is that 
of a spherical emitter that is completely optically thin 
and homogeneous. In this case, the path length of the 
section of the ray that is inside the source is a proxy 
for the brightness received by the corresponding cam- 
era pixel. Figure 4 shows normalized image maps of the 
path length taken through the source, for two spherical 
volumes (emitter proxies) of radii 15 and 5 Rg. Re- 
calling our discussion in Section 2, we see that, for the 
larger emitter, the darkening factor is roughly 1/2, as in 
the Newtonian case; for the smaller emitter, the path- 
lengthening effect dominates, instead. 

Figure 5 shows a similar path-length map, for a vol- 
ume with outer radius 20 Rg and inner radius 10 Ra 
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Figure 4. Upper panels: maps of null-geodesic path length inside a spherical volume of radius 15 Rg (left) and 5 Re (right). The axes 


of the upper panels are measured in Rg. Lower panels: central, horizontal cross-sections of the upper panels; the horizontal axis represents 


the null geodesics’ impact parameter, while the vertical axis represents the normalized path length. Note the step-like dip in intensity that 


occurs at the boundary of the shadow. 


(so that the inner 10 Rg are evacuated). In this case, 
no path-lengthening effect is observed, because the re- 
gion in which that effect occurs (near the boundary of 
the unstable-photon region) is evacuated. Thus, we see 
that in cases like Fig. 5, the blocking effect dominates 
while the path-lengthening effect (and, correspondingly, 
the ‘photon ring’) vanish; in other cases, such as the 
right panel of Fig. 4, it is the path-lengthening effect 
that dominates. We note that this is completely analo- 
gous to the models depicted in Fig. 4 of Narayan et al. 
(2019), which consists of a spherical accretion flow with 
a spherical, inner ‘cutout’ of varying radii; in this model, 
too, the shadow remains constant, even when the cutout 
radius is varied, and the photon ring becomes invisible. 
This provides a counterexample to the claim by Gralla 
(2020) that the photon ring is generic (and that the 
shadow is not). 


Figure 6 shows a path-length map for a volume with 
external radius 20 Rg, but for a Kerr spacetime with 
a = 0.9375. In this case, the non-zero spin gives rise to 
frame-dragging, and the characteristic asymmetric ap- 
pearance of the black-hole shadow. Note, in particular, 
how the path-lengthening effect is suppressed on the im- 
age’s left-hand side. This is because rays that are com- 
ing toward the observer on this side of the black hole 
receive a ‘boost’ due to frame dragging - these rays do 
not plunge, but are quickly ‘swept clear’ of the unstable- 
photon region due to the rotation of spacetime. 


3.2. Toroidal and cylindrical emitters 


In order to understand the effects on the shadow of 
breaking spherical symmetry, as well as the influence 
of the emission-region geometry (i.e., the astrophysical 
circumstances), we now turn to toroidal and cylindrical 
emitters. These models break spherical symmetry be- 
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Figure 5. As Fig. 4, but with an outer volume radius of 20 Rg 
and an inner volume radius of 10 Rg (i.e., no emission originates 
from within the inner 10 Ra). 


cause they carry a non-zero angular momentum, mak- 
ing them more realistic representations of astrophysi- 
cal emission regions. Such emitters are geometrically 
thick and optically thin, corresponding to radiatively 
inefficient, low-density, high temperature emission re- 
gions, which occur around low-luminosity AGN in accre- 
tion flows and jets (Falcke 2001). Because the accretion 
flow carries angular momentum, the effect of relativistic 
boosting becomes apparent for these models, increasing 
the intensity of the side of the torus (or cylinder) on 
which the flow moves toward the observer and diminish- 
ing the intensity of the receding side (the effect vanishes 
when the emitter is seen face-on, and is strongest when 
the observer is in the black hole's equatorial plane). It 
is interesting to note that, despite the absence of the 
path-lengthening effect on that side of the black hole 
(see Fig. 6), this region is still the brightest region in 
an image of a rotating black hole, due to the relativistic 
beaming of the accretion flow. 

In our analytical model, the accretion flow's local ve- 
locity vector is constrained to be tangential to a circular 


-20 -10 0 10 20 


Figure 6. As Fig. 4, but for a Kerr spacetime with a — 0.9375 
and an outer volume radius of 20 Re. 


path. The velocity is then scaled as follows: 
| do | 1 (6) 
E dt ^ SBB Fa 


which, in combination with the demand that the vector’s 
norm equals -1 and a choice of rotation direction (pro- 
grade with respect to the black hole’s sense of rotation 
in our case), determines the vector elements. Only in a 
section of the equatorial plane do the paths implied by 
these tangent vectors correspond to timelike geodesics 
(the stable, circular orbits just outside of the black hole’s 
innermost stable circular orbit (ISCO)). Thus, other re- 
gions of the flow must be affected by pressure gradients 
or magnetic fields. 

The emission coefficient of our toroidal model is pro- 
portional to a Gaussian profile defined as follows: 


(r— = 


jy (r) = exp |— 7 
jv (£) P| E (7) 
where rs is the location of the torus center in the r,0 
plane (since it is in the equatorial plane, 0, = 7/2), and 


o is the parameter that sets the torus’ minor radius. 


Using this model, we now plot several intensity maps 
for qualitatively different scenarios, to illustrate how 
the different physical processes behind the black-hole 
shadow contribute. 

Figure 7 shows a torus with major radius 16 Rg and 
minor radius 3 Rg, viewed from the equatorial plane. 
In this model, the emission region lies well outside the 
black hole’s ISCO, and demonstrates the ‘evacuation’ 
effect discussed in Bronzwaer et al. (2020), causing an 
enlarged CBD around the black hole’s true shadow, al- 
though the sharp-edged optical signature at the shadow 
remains. Figure 8 shows the effect of revolving this emis- 
sion region around the image X-axis, i.e., summing the 
images from all inclination angles; this effectively ‘fills 
in’ the evacuated regions, but the step-like dip in bright- 
ness at the shadow remains constant, thus producing a 
clearly observable shadow. 


-20 -10 0 10 20 


Figure 7. Intensity map of our toroidal emission model with 
Rmaj = 16 Rg and Rmin = 3 Ra, in a Schwarzschild spacetime, 
for a distant observer situated in the black hole’s equatorial plane. 


Figure 9 shows a Kerr black hole (a = 0.9375) sur- 
rounded by a toroidal flow with Rmin = 1 Re on an orbit 
that lies 1 Rg outside of the ISCO. Note that the high 


Figure 8. The model shown in Fig. 7, revolved around the image 
X-axis (a summation of intensity maps for all inclination angles). 


-20 -10 0 10 20 


Figure 9. Intensity map of our toroidal emission model with 
Rmaj = Risco + 1 Re and Rmin = 1 Ra, in a Kerr spacetime 
with a = 0.9375, for a distant observer situated in the black hole’s 
equatorial plane. 
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spin causes a strong relativistic boosting effect and thus 
a small emitter size, as we see only the boosted region. 
Compared to the Schwarzschild case, the ISCO lies much 
closer to the unstable-photon region for a high-spin black 
hole, and thus, the evacuation effect is reduced, and a 
clear black-hole shadow is observed for co-rotating disks 
around high-spin black holes. Interestingly, Figs. 7 and 9 
are rather similar in appearance to the equivalent SANE 
disk models discussed in Bronzwaer et al. (2020). 

In order to explore the range of possible models be- 
tween the geometrically thin and geometrically thick 
regimes, Fig. 10 shows the transition from a thin to a 
thick torus, illustrating the decreasing severity of the 
evacuation effect for thinner torii as the disk gets thicker. 
For very thick disks, a clear shadow can be seen, rem- 
iniscent of the spherical case, as the blocking effect is 
on full display. Thinner emission regions, on the other 
hand, contain relatively more evacuated regions, inter- 
fering with the blocking effect and rendering the shadow 
less easily observable; a step-like dip in brightness is still 
observed, but the CBD has a more complicated shape. 

From the preceding examples, we see that in the case 
of toroidal emitters, the visual signature (step-like dip in 
brightness) that is coincident with the BHS arises due to 
a complex interplay of gravitational lensing (causing its 
exaggerated size), the path-lengthening effect (increase 
of the brightness along its edge), and the blocking effect 
(relative darkening of rays that intersect the unstable- 
photon region). It is remarkable that these effects all 
contribute to a visible signature that lies on the edge 
of the shadow (as seen by a distant observer). For a 
more detailed look at GRMHD-based toroidal emission 
models, see (Bronzwaer et al. 2020). 

Figure 11 shows our cylindrical emission model, which 
is a simple model for a disk-jet accretion flow. The cylin- 
der has a hollow core; its radius is 9 Re and its wall 
thickness is 6 Rg. As in the toroidal case, the cylin- 
der is given a purely azimuthal velocity vector which is 
scaled by Eq. 6. Lensed images of the emission region’s 
hollow inner core can be observed just above and below 
the shadow. Note, however, that the step-like dip in 
intensity at the BHS is visible in each case; this model 
contains only tiny and well separated evacuated regions, 
compared to the toroidal models, so that essentially a 
complete shadow is observed. 


3.3. Thin-disk emitters 


Emission regions with a higher density and radiative 
efficiency and a lower temperature than the RIAF flows 
discussed in the preceding section may form in, e.g., 
black hole binaries and QSO’s. A model was devel- 
oped to describe such flows (Shakura & Sunyaev 1973; 


Novikov & Thorne 1973); it consists of a geometrically 
thin, optically thick disk confined to the equatorial plane 
and terminating abruptly at the ISCO. The velocity vec- 
tor is identical to the previous section, so that disk el- 
ements move along circular, timelike geodesics. As the 
disk is geometrically thin and optically thick, its associ- 
ated emission coefficient is a delta function along a ray, 
and the ray terminates upon encountering the disk. 

Figure 12 shows intensity maps of this model for an 
observer-inclination angle of 75deg and two different 
black-hole spins, as well as a Newtonian analog of the 
model, which is included for comparison; note the dif- 
ferent length scale of the Newtonian image, which indi- 
cates the difference in apparent size of the source due 
to the absence of gravitational lensing. Regarding the 
BHS for these models, as in the case of toroidal emit- 
ters, the evacuation and obscuration effects are apparent 
(the evacuation effect is rather prominent for this model, 
as we exclude any emission from within the ISCO, a 
crude assumption which makes it a sort of worst-case 
scenario). 

In this model, the blocking effect and the path- 
lengthening effect are entirely absent, as there is no ge- 
ometrically thick, optically thin flow through which to 
measure path length in the first place. Thus, no step-like 
drop in brightness occurs at the shadow, and the appear- 
ance of the CBD is much more model-dependent than in 
the case of an optically thin, geometrically thick emis- 
sion region. However, the gravitational lens effect still 
displaces the emission coming from behind the unstable- 
photon region (with respect to the observer) so that it 
appears to come from outside of the unstable-photon 
region. Thus, as the disk’s inner radius approaches 
the unstable-photon region, a CBD may be observed 
that is nearly coincident with the BHS. This occurs at 
high spins, for which the ISCO and the unstable-photon 
region both shrink down, approaching the size of the 
event horizon for a — 1 (Bardeen et al. 1972). At 
low spins, or particularly retrograde spins, the ISCO is 
pushed further out, and the evacuation effect described 
in Bronzwaer et al. (2020) exaggerates the size of the 
observed CBD. The obscuration effect described in that 
same paper also occurs in the thin-disk model, for most 
inclinations. However, the evacuated region differs from 
the true shadow’s interior, because the celestial sphere is 
visible in the former, and, depending on the observing 
frequencies, astrophysical objects (e.g., distant AGN) 
will be visible there, thus rendering the two regions dis- 
tinguishable, at least in theory. Of course, in practice, 
the brightness of the lensed objects that appear in the 
evacuated region is much lower than that of the accre- 
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Figure 10. Intensity maps for a toroidal emitter with Rmaj = 9 and Rmin = 3, with a z-scaling factor of 0.5 (left panel), 1 (center panel), 


and 2 (right panel). 


Figure 11. Intensity map of our cylindrical emission model with 
a radius of 9 Rg and a wall thickness of 6 Rg (so that the inner 
core is evacuated). 


tion disk, and it may also be difficult to disentangle this 
emission from any foreground emission. 


4. DISCUSSION AND CONCLUSIONS 


A black hole’s shadow is the projection of a its 
unstable-photon region on an observer’s sky. If the black 
hole is surrounded by an optically thin, geometrically 
thick emission region, a step-like dip in brightness will 
be observed that is coincident with the shadow. This 
dip may or may not be coincident with the overall CBD 
that is observed (which may itself have a complicated 
structure, see, e.g., Chael et al. (2021)). The purpose of 
this work was to get a more intuitive understanding of 
the shadow as an astrophysical observable. 


In the case of spherical emission regions, the primary 
mechanism behind the black hole shadow is the blocking 
effect; the shortening of path lengths for null geodesics 
that intersect the photon sphere and terminate at the 
event horizon, versus null geodesics that do not inter- 
sect the photon sphere. For these spherical models, the 
shadow may be readily identified by an edge-detection 
scheme (Psaltis et al. 2015) or some form of crescent 
shape fitting (Kamruddin & Dexter 2013; Event Hori- 
zon Telescope Collaboration et al. 2019b). 

Besides the blocking effect, the path-lengthening ef- 
fect occurs, whose observational signature is a relative 
enhancement of brightness at the edge of the black-hole 
shadow. We note that this effect occurs both inside 
and outside of the unstable-photon region (and thus 
the shadow); the cusp-like increase in brightness caused 
by the path-lengthening effect peaks on the shadow’s 
boundary, and decays in both the outward and inward 
directions. However, when the observable effects of the 
blocking effect and the path-lengthening effects are com- 
bined, the part of the cusp that is interior to the black- 
hole shadow appears reduced in brightness (see the pro- 
file of the left panel of Fig. 4), so that the rays that 
travel exterior to the unstable-photon region dominate 
the image, and the shadow’s interior is overall dimin- 
ished in brightness. These strongly bent light rays, 
which travel near the edge of the unstable-photon re- 
gion (but which are outside of it), produce a complex 
system of lensed images of the black hole’s environment. 
This is often called the photon or light ring, and it has 
been invoked for tests of gravity (Luminet 1979; Fal- 
cke et al. 2000b; Johannsen & Psaltis 2010a; Johnson 
et al. 2020; Gralla et al. 2020). The photon ring con- 
tains multiple projections of the surrounding emission 
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Figure 12. Intensity maps of the thin-disk model by Novikov & Thorne (1973), for an observer-inclination angle of 75 deg, and a black- 
hole spin of a = 0 (left panel), and a = 0.9375 (center panel). Right panel: Newtonian analog of thin-disk model, also viewed at i = 75 deg. 


region, and approaches the shadow asymptotically the 
more often a light ray goes around the black hole (thus 
approaching the peak of the cusp); Moscibrodzka et al. 
(2016) show that, in the case of M87*, the lensed, first- 
order image of the counter-jet can actually dominate 
the image of the source (see also Appendix A of Dav- 
elaar et al. (2019)). We have shown that a step-like 
dip in brightness at the shadow, and even a CBD that 
is coincident with the shadow, can arise independently 
from the path-lengthening effect (and thus the photon 
ring). This will happen, for example, if a volume out- 
side the black hole’s ISCO is evacuated, as in Fig. 4 of 
Narayan et al. (2019) and Fig. 5 of this work. In certain 
models, the CBD can be ‘multi-tiered’; a step-like dip 
may occur at the shadow, along with an ‘inner shadow’, 
which is a smaller, smooth-edged dip that occurs at the 
projection on the observer’s sky of the intersection of 
the black hole’s event horizon with its equatorial plane 
(Chael et al. 2021). 

Neither the blocking effect nor the path-lengthening 
effect is relevant with regard to the appearance of geo- 
metrically thin emission disks. The dark region (CBD) 
inside of a thin disk can be larger or smaller than the 
BHS, as one observes a combination of the shadow and 
the evacuated region in the emission disk. This ambigu- 
ity disappears when the emission regions transits from 
a planar to a spherical structure (See Fig. 10). 

We note that the appearance of a sharp-edged CBD 
that is coincident with the projection of a dark object 
(such as a black hole’s unstable-photon region) is itself 
not unique to GR. A perfectly light-absorbing object 
in flat (Newtonian) spacetime and surrounded by an 
optically thin emission region would produce a similar 
CBD. The everyday example we gave is that of a kettle 


in a campfire, or of a wick in a candle. However, the 
shadow in GR is exaggerated in apparent size by a fac- 
tor of about 2.5 due to gravitational lensing. Moreover, 
the image can exhibit a bright photon ring due to path 
lengthening in the emission region, if the emission re- 
gion (nearly) intersects the photon sphere. Both effects 
are due to light bending. The darkening itself, how- 
ever, implies either the presence of an event horizon - a 
perfectly absorbing surface with zero reflectivity at all 
wavelengths - or of another type of surface which must 
nevertheless lie close to the event horizon, and thus is 
hidden from outside observers due to the strong grav- 
itational redshift. This is in contrast to the everyday 
Newtonian example of the kettle in a campfire, where 
the initially dark kettle would heat up and eventually 
be seen in emission. The absence of such thermal emis- 
sion has been used to infer the presence of a horizon 
in stellar-mass black holes (Narayan et al. 1997) and 
in M87* (Event Horizon Telescope Collaboration et al. 
2019b). 

In observations, the darkness of the observed CBD 
is limited by the dynamic range and resolution of the 
observing instrument, as well as the accuracy with which 
foreground emission can be subtracted. Therefore, the 
presence or absence of a black hole's event horizon can 
only be verified to within certain limits. However, since 
the emission region is close to the black hole, any light- 
emitting matter moves at a significant fraction of the 
speed of light, whether falling into the event horizon, 
escaping in a plasma jet or wind, or rotating on bound 
orbits. The dynamical time scale for these processes 
range from milliseconds for stellar-mass black holes to a 
couple of days for a supermassive black hole like M87*. 
Hence, integrated and averaged over human time scales, 


any irregularities in the foreground emission should be 
smoothed out. 

In case the black hole is not completely surrounded 
by an emission region, an evacuated region (Bronzwaer 
et al. 2020) may be observed around the shadow. How- 
ever, unlike the shadow itself, the darkness in these 
patches is of a fundamentally different nature. It will 
not be completely dark and smooth, but will contain a 
static lensed image of the universe. The same applies for 
the criticism of Bardeen (1974) by Gralla et al. (2019), 
who states that a black hole in front of a uniformly radi- 
ating surface, e.g. a stellar disk, appears larger than the 
projection on the observer’s sky of the unstable-photon 
region. However, in this case, too, the ‘dark’ region sur- 
rounding the true shadow is of a fundamentally different 
nature than the shadow itself, as it, too, is filled with 
emission from lensed images of distant sources in the 
entire universe. Unlike any rapidly moving foreground 
emission, this background image will remain static, al- 
though it will also be extremely faint. 

Some compact objects without horizons (but with 
highly redshifted surfaces deep within their photon 
spheres) can also show shadows. However, such objects 
may be distinguishable from black holes through either 
re-radiation (Broderick & Narayan 2006) or reflection 
of emission. For example, in the absence of a photon- 
capture cross section, the central dark region in the bo- 
son star case (Vincent et al. 2016) is simply a lensed 
image of the star’s central low-density region (Olivares 
et al. 2020). Light could pass through the boson star it- 
self. In rotating boson stars, plunging geodesics can be- 
come chaotic, and are reflected outwards towards infin- 
ity (Cunha et al. 2016; Olivares et al. 2020). One would 
expect to see a highly distorted, lensed image of the 
universe in lieu of a perfectly dark shadow, for such ob- 
jects. Naked singularities could, in principle, also show 
a shadow (Shaikh et al. 2019) as they are matched to a 
Schwarzschild solution on their exterior, and can possess 
a photon sphere. 

These are just some examples of possible ambigui- 
ties that will generally exist for any non-horizon object 
whose radial extent makes it fall within the ‘twilight 
zone’ between the true horizon and the photon sphere. 
One can always think of a deviation small enough to 
escape detection. In any case, black holes are the most 
straightforward and general astrophysical interpretation 
for a range of phenomena. For the case of supermas- 
sive black holes, the shadow was properly predicted well 
before the observations with up-to-date computer sim- 
ulations (Falcke et al. 2000a; Broderick & Loeb 2006; 
Dexter et al. 2012; Moscibrodzka et al. 2016) and then 
observed. We have explained here that indeed a black 
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hole’s shadow is a clearly defined and robust feature as 
long as the emission region surrounding the black hole is 
geometrically thick and optically thin. With these basic 
assumptions, the exact intensity distribution or details 
of the astrophysics and particle heating are largely ir- 
relevant. Indeed, for jet sources like M87* it is well 
established that the core becomes more optically thin 
at higher frequencies, and that the emission is expected 
to be spatially extended in all directions. The measured 
constancy of the shadow diameter of M87* over a decade 
Wielgus et al. (2020) also renders an interpretation of 
the image as incidental, transient astrophysical circum- 
stances of the accretion disk/jet system unlikely. 

History has shown that it is impossible to derive 
physics from astrophysics without understanding astro- 
physics. The crucial question is not whether current or 
future EHT results depend on astrophysics, but rather 
how significant the uncertainty due to unknown astro- 
physics is. The EHT has gone to great lengths to quan- 
tify this uncertainty using a wide range of self-consistent 
numerical simulations that take the behaviour of plasma 
radiation into account. We have come a long way to 
understanding how black holes accrete and how they 
launch jets. Simple toy models, therefore, can only be 
used to guide our intuition and help us to explain what 
we are seeing. Nonetheless, care has to be taken that 
these toy model at least roughly reflect the actual astro- 
physical situation, which is not entirely arbitrary any- 
more. 

In summary, we reiterate that the black-hole shadow 
is a robust and significant observational feature, which 
speaks of the underlying metric and the darkness of 
the event horizon, and which is independent of the 
phenomenon of photon rings. The first detection of a 
shadow in M87* confirms a basic prediction made for 
low-power black holes. Of course, this is only a starting 
point: better tests will require more theoretical investi- 
gations, more observations, and better experiments, in- 
cluding space experiments. Uncertainties due to source 
variability will decrease by repeated EHT observations, 
and will sharpen constraints on key parameters (see, 
e.g., Roelofs et al. 2021, Psaltis et al. 2020). Given con- 
tinuous rubbing of the coin, a more solid outline of the 
underlying metric will appear, and we will be able to 
read the mint markings of black holes more clearly in 
the future. 
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